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The dynamical response of a trapped Bose-Einstein condensate (BEG) is formulated 
consistently with quantum field theory and is numerically evaluated. We regard the BEG as 
a manifestation of the breaking of the global phase symmetry. Then, the Goldstone theorem 
implies the existence of a zero energy excitation mode (the zero-mode). We calculate the 
effect of the zero-mode on the response frequency and show that the contribution of the 
zero-mode to the first excitation mode is not so important in the parameter set realized 
in the existing experiment. This is the reason that experimental results can be described 
using the Bogoliubov prescription, although it breaks the consistency of the description in 
quantum field theory. 

§1. Introduction 

Since the experimental realization of the Bose-Einstein condensate (BEC) in 
trapped dilute atomic g large number of works associated with trapped BECs 

have been reported. The majority of these works are devoted to the analysis of the 
Gross-Pitaevskii (GP) equation.^) The GP equation has been successful in account- 
ing for experimental results at very low temperature. However, the GP equation 
is not capable of describing the BECs as the temperature is raised, because, for 
instance, we should take the effect of the non-condensate fractions into account. 
Moreover, we can continue to fine-tune experiments on trapped BECs, changing the 
scattering length between atoms, by utilizing Feshbach resonance.^) Thus, trapped 
BECs afford an opportunity for testing quantum field theory (QFT), which is still 
incomplete. To this end, we need to analyze trapped BECs beyond the GP equation. 

There are several theoretical works in which the behavior of BECs at finite tem- 
perature both with^)"^^) and without a trapping potential^^^"^^-' is studied. In par- 
ticular, the results given in Ref. 10), as well as Refs. 11)-13), are in good agreement 
with the experimental value^^^ of the excitation frequencies from the low tempera- 
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ture region to the region near the transition temperature. The authors of Ref. 10) 
use Zaremba-Nikuni-GrifHn formalism/^) in which the time-dependent GP equation 
couples with the scmiclassical Boltzmann equation. In Refs. 11)-13), the authors 
apply a time-dependent extension of the perturbation theory. 

The quantization scheme of a trapped BEC is, however, an open question. From 
the viewpoint of QFT, the BEC is a manifestation of the spontaneous breakdown 
of a global phase symmetry. In this picture, the Goldstonc theorem requires the 
existence of a mode with zero excitation energy (the zero-mode) for consistency with 
the fundamental concepts of QFT. However, a consistent quantization scheme is not 
easy to construct, because of, for instance, the problem of the infrared divergence. 
Usually, when the Bogoliubov prescription is employed, the zero-mode is ingored, 
but the naive application of the Bogoliubov prescription gives rise to inconsistency, 
because the precise canonical commutation relations cannot be satisfied. Thus, we 
cannot ignore the zero-mode in testing the predictions of quantum field theory in 
trapped BEC experiments. 

Recently, a consistent quantization scheme was proposed by two of the present 
authors.^^^ They introduced an infinitesimal breaking term of the global phase sym- 
metry. Here, the breaking term plays roles as a regulator of infrared divergences 
and enables us to construct the generalized Bogoliubov transformation (GBT) in- 
cluding the zero-mode. It has been proven that this mode corresponds to the 
Nambu-Goldstone (NG) mode by demonstrating that the Ward-Takahashi (WT) re- 
lations^'^)' ^-"^^ are satisfied. Realization of unitarily inequivalent vacua that belong 
to order parameters with different phases was shown in Ref. 22). The existence of 
such vacua implies that the mechanism of the the spontaneous symmetry break- 
down works even in traps. There is another consistent field theoretical approach, 
in which the operators corresponding to the zero-mode are given by "quantum co- 
ordinates". ^^^'^^^ The relation between the zero- mode approach and the quantum 
coordinate approach has been revealed, but it has not been confirmed that the 
WT relations are satisfied in the latter approach. It should also be mentioned that no 
definite criterion has been found for choosing quantum states acted on by quantum 
coordinates. Effects of the zero-mode, as well as modes represented by the quantum 
coordinates on the condensate density have been evaluated, and, unfortunately, 
they were found to be too small to be detected in present-day experiments. 

The aim of the consistent canonical studies cited above is to elucidate the equi- 
librium properties of trapped BECs.^'^-*"^^) In this paper, we apply the consistent 
quantization scheme and study the effect of the zero-mode on the dynamical response 
of a trapped BEC at finite temperature. 

This paper is organized as follows. In §2, we explain the quantization of a trapped 
BEC including the zero-mode. In §3, the dynamical response of the BEC under an 
external field is derived, and the results of numerical calculations are reported in §4. 
Section 5 contains concluding remarks. We give a brief summary of the mathematics 
regarding the GBT including the zero-mode in Appendix A. 
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^2. Quantization scheme consistent with the Goldstone theorem 



In this section, we give a brief review of the quantization of a static trapped 
BEC including the zero-mode, following Refs. 21), 25) and 26). 

We adopt the following action, which describes the dynamics of a neutral gas 
system trapped in a potential V: 



dtd^x 



^Hx){T - K -V + n)^{x) -^^^'^{x)ij'^{x) . (2-1) 



Here, T = i-^, K = — 2^V^, and /x and g denote the chemical potential and the 
coupling constant, respectively. We choose the harmonic trapping potential as ^ = 
T Y^i=i ^i^h where is the trapping frequency along the Xi direction. The symbol 
X represents {x,t). For simplicity, we set ^ = 1 throughout this paper. 
This action is invariant under the global phase transformation 



tp{x) 



ip{x) 



and 



-IT] 



ip^{x) 



(2-2) 



where rj is an arbitrary constant phase. The Bose-Einstein transition can be regarded 
as the spontaneous symmetry breaking associated with this phase transformation. 

Usually, the action (|2-ip is the basis of the quantization of trapped BEC systems. 
However, in analogy to the Ising model, we add the following artificial symmetry 
breaking term to Eq. ()2-l() : 



Z\5[^,^t] dtd^x 



%{x)il;{x) + e'^v{x)ip'^{x) 



(2-3) 



Here, e denotes a typical energy scale of the system determined by trapping frequen- 
cies. The infinitesimal parameter e is taken to be vanishing at end of the calculation. 
This term is necessary to control the infrared diver gence.^^^ The order parameter 
e^^v{x) is determined by a consistent procedure. 

Now we consider canonical formalism. The quantized field satisfies the CCRs 



^{x,t),^j\x'A =<)-(3)( 



[X — X 



^{x,t),^{x',t) \ = \ij\x,t),^\x',t) 



0, 



(2-4) 
(2-5) 



where the symbol " indicates that quantities with it are regarded as operators. Cor- 
responding to the breakdown of the global phase symmetry, we divide the original 
field il^{x) into the classical and quantum parts, as 



^Ij{x) = e'"v{x) + e^V(a;) • 



(2-6) 



Here, the time-independent real function v{x) is a condensate field that breaks the 
global phase symmetry defined by Eq. H2-2|) . Thus, v{x) is the order parameter of 
the BEC transition. It should be noted that |f(a;)p gives the number density of 
the condensate particles. It is assumed that the real number 9 is independent of 
both time and space, corresponding to the situation without vortices. In this case, 

*-* When we consider a vortex, the phase depends on the space coordinates. 
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without loss of generality, we can set 9 = 0. From the CCRs of il^ and ifj^ , (|2-4j) and 
H2-5|l . the CCRs of tp and (p^ are given by 



[(p{x,t),(p{x',t)] = ip\x,t),(p^{x',t) 



0. 



Calculating the Hamiltonian corresponding to the total action S 
substituting Eq. H2-6p into it, we obtain 



where 



Ho 
Hi 



d^x 



H = Hq + Hi , 



ip\K + V - n)ip + ^-{i^^ip + ip'^ + 



t2^ 



(2-7) 
(2-8) 

AS, and 
(2-9) 

(2-10) 



d^x 



v{K + V — ^ + gv^ — ee)(p + tp^{K + V — + gv^ — ee)v 



+ gv{ip^Lp^(p + p^p0) + ^^^0^00 



(2-11) 



Here, we have dropped the constant term independent of and 0^^ . 

Now we adopt the Hartree-Fock-Bogoliubov approximation, in which the third- 
and fourth-order terms of and 0^^ are replaced as follows: 



0'^0'^0 > {0'^0'^)0 + 20'^{0'^0), 

0'^00 > 2{0^0)0 + 0^00), 

0'^0'^00 > ^{0'^0)0'^0 + {0^'^)0'^ + {0'^)0^'^. 



(2-12a) 
(2-12b) 
(2-12c) 



Here, ( ) represents the thermal expectation value, which will be determined self- 
consistently using the mean- field approximated Hamiltonian Hmfa, and which is 
defined through the relation (O) = Trle'f^^^^^dj/Z with Z = Tr[e"^''^MFA]^ where 
O is any operator. We further employ the Popov approximation and drop the terms 
proportional to {0"^) and (c^^^)^^) in order to ensure a gapless spectrum.*^ Finally, 
the Hamiltonian is reduced to 



Hmfa = j d^x ^0^K + V - fi + 2gv^ + 2g{0^0))0 + ^ 

+ v{K + V - n + gv^ - ee + 2g{0^0))0 
+ 0\K + V - ^ + gv'^ - ee + 2g{0^0))v 



9 -t2 2 



(2-13) 



*' For a homogeneous system, it has been proven that in the thermodynamic Hmit, if one 
performs a second-order perturbative calculation, the Hugenholtz-Pines (HP) theorem, which 
ensures the global gauge invariance, holds. There is yet no similar proof for a trapped system. 
Recently, it has been shown'^^' that the HP theorem holds for the loop expansion in the case of 
a finite volume system without a trapping potential but with periodic boundary conditions even 
without taking the thermodynamic limit. 
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The terms first-order in (p and ip'^ siiould vanish. Thus, the condensate v satisfies 

(K + V - fi + gv'^ - 66 + 2g{ip^p))v = 0. (2-14) 

This equation is identical to the GP equation^^ in the limit e ^ 0. 

Because the von Neumann theorem does not hold in QFT, there are many Fock 
spaces that are orthogonal to each other,^^) and we need to choose a physical Fock 
space among them. Here, we adopt the procedure proposed by two of the present 
authors, in which the physical Fock space is defined by the vacuum of annihilation 
operators defined by diagonalizing i^MFA using the GBT.^^^'^^^'^^) 

First, we prepare an orthonormal set {wn{x)} obtained by solving the differential 
equation 

(^K + V-n + gv'^ + 2g{(p^ip)^ Wn{x) = (e„ + ee)wnix). (2-15) 
The orthonormal and completeness conditions of this set are 

j d^XWn{x)Wn'{x) = 5nn', (2-16a) 

oo 

^Wn{x)wn{x') = 5^^\x - x') . (2-16b) 

n=0 

Then, the field operators (p{x) and (p\x) are expanded as 

oo 

"^l^) = ^^nWCn{x) - h\wsn{x)] , (2-17a) 
n=0 

oo 

^\x) = ^[hiwcn{x) - hnWSn{x)] , (2-17b) 



n=0 



where 



WCn{x) = ^ CnmWm{x), (2-18) 

m=0 

oo 

WSn{x) = ^ SnmWm{x). (2-19) 
m=0 

The matrices C and S are introduced in association with the GET that diagonalizes 
Hmfa- (See Appendix A for details.) Note that C and S are singular in the limit 

e and give rise to an infrared divergence. 

Finally, the Hamiltonian Hmfa is diagonalized as follows: 



oo 



^MFA = ^nbik + const , (2-20) 

n=0 

where the quasiparticle energies depend on {wn{x)}, v{x) and {e„}, and its 
expression is given in Appendix A. It should be noted that the energies £'„ are finite 
even after taking the limit of e — >^ 0. Finally, the physical vacuum is defined by 

S„|0) = 0. (2-21) 
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The creation and annihilation operators satisfy the bosonic commutation rela- 
tions 

[bn,bl]=5nn', (2-22) 

with all other commutators vanishing. Thus, the quantization scheme explained 
here satisfies the CCRs (|2-7() and ()2-8|) precisely, in contrast to the Bogoliubov pre- 
scription. This is because the mode associated with Eq is not included in the field 
expansion in the Bogoliubov prescription. 



§3. Dynamical treatment of a system under an external perturbation 

To this point, we have studied the quantization scheme in a static system. In 
this section, we discuss the dynamical behavior of a EEC. 

To investigate the response to an external field, we introduce the time-dependent 
external perturbation 



3 

m 



H^^{t)= j d^xij\x)5V{x,t)i^{x), (3-1) 

where 

5V{x,t)='^Y.^u:Ut)xl (3-2) 

i=l 

Then, the total trapping potential is given by 

VUx,t) = V{x) + 5V{x,t). (3-3) 
The Heisenberg equation of motion with a time-dependent perturbation is given by 

^^'(/'ex = {K + Vc^{x, t) - /i - ee + gV'exV'ex)'0ex, (3-4) 

where ip^^ denotes the boson operator subject to the perturbation. As in the pre- 
ceding section, we separate ipe-^ into two terms as follows: 

i)e^{x) = Vex{x,t) +ip{x). (3-5) 

We further employ the Hartree-Fock-Bogoliubov-Popov (HFBP) approximation. Then 
the GP equation becomes 

d 

^^^^ex = (i^ + Vex - Ai - ee + ffbexP + 2g{(p'^Lp)cx)Vex, (3-6) 

where |wexP is the number distribution of the condensate particles. Now, the order 
parameter has an explicit time dependence, in contrast to the static case. 

However, Eq. (|3-6|) is still difficult to solve, because we must redefine a physical 
Fock space with the time evolution of the condensate. If we are interested in low 
temperature and/or weak coupling systems in which fluctuations are expected to be 
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small, however, we can further linearize the equation about the equilibrium values. 
Thus, the condensate field and the number density of the excited particles, {ip^ip)ex, 
are expanded about the static quantities defined in the preceding section, 

Vex(x) = v{x) + 6v{x), (3-7) 
{^^ip),x = {ip^^) + 5{^^^). (3-8) 

Keeping the first-order terms 6V{x), 6v{x) and 6{(p^(p) and ignoring the higher-order 
terms, we obtain the linearized equation. It should be noted that the linear approx- 
imation with the HFBP approximation can be regarded as the time-dependent lin- 
earized Hartree approximation, and this approximation is known to be equivalent to 
the random phase approximation (RPA). Thus, this approach contains fluctuations 
beyond the mean-field approximation. To obtain further simplicity, we may ignore 
the term proportional to 6{ip^ip), because the quantity {(p^(p) is already very small. 
Finally, we obtain 

i—5v{x) = {k + V - 11 + 2gv'^ + 2c,((^V)) 6v{x) + gv^6v*ix) + 6V{x)v . (3-9) 

Note that the phase of the condensate field is time dependent. Thus, this equation 
can in principle describe vortex formation. 

The effect of the zero- mode is included in Thus, it possesses infrared 

divergence, because of the sing ularity of S and C'.2i)'25),26) ^i^jg difficulty, 

we employ the renormalization criterion according to which all terms proportional 
to u'^{x) are ignored. It should be noted that we still keep the terms proportional 
to ui{x)uq{x) {I 7^ 0) in the limit of e — > 0. Thus, the nontrivial effect attributed to 
the zero- mode is still taken into account. 

§4. Numerical study 

In this section, we report the results of our numerical study of Eq. (|3-9|) . In this 
study, we chose the parameter values so as to reproduce the experimental situation.^'*) 

4.1. Thermal particle density 

According to Eq. ()3-9|) . the time evolution of the condensate is affected by the 
thermal non-condensate particle density {(p'^ip). First, we study the effect of the 
zero-mode on the non-condensate particle. For simplicity, we treat a spatially one- 
dimensional system and consider the BEC of rubidium atoms, ^''Rb, whose mass 
is m = 1.42 X 10~^^kg, and set the frequency of the trapping potential as uJi = 
200 X 27rHz. Here, we set the parameter values as T/T^ = 7.3 x 10~^, g = 0.02 in 
oscillator units (HO units), with the condensed atom number A'c = 1000. Here, Tc is 
the transition temperature of an ideal Bose gas in a one-dimensional system, whose 
value is T^ = = 8.3 x 10~^K, according to the estimation given in Ref. 34), 

where /cb is the Boltzmann constant and is the total number of atoms, which in 
our calculations is A = 1068. The numerical calculations were carried out in the 
following steps. First, we diagonalized the mean field Hamiltonian H\i\-ph (|2-13() under 
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Fig. 1. Thermal particle densities of the one-dimensional system, with the parameter values T/Tc = 
7.3 X IQ-^ g = 0.02 (HO units). The solid curve represents the result with the zero-mode, while 
the dotted curve represents the result without it. 

the condition ()2-14|l to define the physical Fock space. Accordingly, we estimated 



Figure^plots this thermal particle density, calculated both with and without the 
zero-mode. The solid and dashed curves represent the thermal particle density with 
and without the zero-mode, respectively. One can clearly see that the density of the 
thermal particles increases around x = 0, because of the zero-mode. This difference 
is, however, very small compared to the condensate number density, because the 
thermal particle density is about two orders smaller than the condensate particle 
density in our calculation. Thus, it is still not clear whether the zero-mode can 
significantly affect the dynamical response of the condensate field. 

4.2. Time evolution of the condensate 

The actual experiment was carried out in a cylindrically symmetric system. 
However, we assume spherical symmetry here for calculational simplicity. Then, 
the static trap potential with frequency to = tOi {i = 1,2,3) and the corresponding 
external perturbation are given by 



the expectation value {^f^'p)- 




respectively. In conformity with the experiment,^^) we employed the following time- 
dependent frequency: 




The quantities A and i? are the driving amplitude and the driving frequency, respec- 
tively. 
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perturb off 

1.775 I ^ ^ ^ n 




1.745 I • • • ^ • • • • ' 1 

5 10 15 20 25 30 35 40 45 50 

time (HO units) 

Fig. 2. Time evolution of the condensate with the parameter values T/Tc = 0.11, A = 0.01, f2 = 1.3 
(HO units) . The expression "perturb off" refers to the time at which the perturbation is turned 
off. After that time, the condensate oscillates sinusoidally with an amplitude of 0.3 % with 
respect to the initial size. 

The numerical calculation was carried out in the following steps. First, we es- 
timated the expectation value {^^0), using the procedure described in the previous 
subsection, and solved Eq. (|3-9() . The parameter values were chosen so as to repro- 
duce the trapped atoms of ^"^Rb, where m = 1.42 x 10~^^kg, to = 200 x 27rHz and 
g = 0.03 in HO units and condensed atom number is Nc = 1000. The temperature T, 
in units of Tc, the driving amplitude A, and the driving frequency f2 are the control 
parameters in the experiment, where Tc = 0.94^A^^/'^ is the critical temperature 

of the ideal Bose gas confined in a three-dimensional spherically symmetric trap.^-* 
The value of Tc is approximately 100 nK in this calculation. 

In Fig. El we plot the time evolution of the size of the condensate, defined by 

The temperature, the driving amplitude and the driving frequency are given by 
T/Tc = 0.11, A = 0.01 and 17 = 1.3 (HO units), respectively. Under the external 
perturbation, we observe beat-like behavior. The external perturbation was turned 
off at about 15.4 ms, which is a typical value in the corresponding experiment.^^^ 
Then, the frequency of the oscillation of the condensate is 2.014 (HO units). This 
frequency characterizes the dynamical behavior of the condensate. We plot the 
numerical result with a different driving frequency, = 2.0 (HO units), in Fig. |21 
After turning off the external perturbation, the frequency of the oscillation of the 
condensate is again given by 2.014 (HO units). It should be noted that the resonance 
gives rise to an increase of the amplitude, and we verified that the amplitude is 
maximized at 17 = 2.014 (HO units). 

The first excitation energy of the quasiparticle is Ei = 2.014 (HO units) [see Eq. 
H2-20() ]. It is thus seen that the frequency of the oscillation is characterized by the 
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Fig. 3. Time evolution of the condensate with the parameter values T/Tc = 0.11, A = 0.01, SI = 2.0 
(HO units). The amplitude of the oscillation is 4.5 % with respect to the initial size. 



first excitation energy of tfie quasiparticle. This result is consistent with that of the 
linear density response. ^^-^ 

4.3. Temperature dependence of the response frequency and the effect of the zero- 
mode 

In the preceding section, we showed that it is possible to extract information 
concerning the first excitation energy of the quasiparticle through the frequency of 
the oscillation. In this section, we investigate the effect of the zero-mode on the 
frequency. 

In Table ^ the temperature dependence of the frequency is described for two 
cases, those in which the frequencies are calculated with the zero-mode and those 
without the zero-mode. The driving amplitude and the driving frequency are given 
hy A = 0.01 and fi = 2.0 (HO units), respectively. The table lists the temperature, 
the frequency with the zero-mode, that without the zero-mode, and the difference 
between these two frequencies. It can be seen that the effect of the zero-mode on the 
response frequency is very small. This is consistent with the fact that the properties 
of a trapped BEG at low temperature are described in the Bogoliubov prescription, 
although it violates the OCRs. It is also seen that the frequency decreases as the 
temperature is increased. It follows that the first excitation energy is lowered by the 
influence of the medium effect. Thus, the energy shift due to the zero- mode becomes 
important as the temperature is increased. However, the magnitude of this shift is 
still very small, and the effect of the zero-mode is negligible for this set of parameter 
values. 

4.4. Coupling dependence of the first excitation energy and the effect of the zero- 
mode 

In the preceding subsection, we showed that the effect of the zero-mode in the 
experiment reported in Ref. 14) is quite small and does not affect the properties of 
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Table I. Temperature dependence of response frequency, calculated with and without the zero- 
mode. 





response frequency 






with zero-mode without zero-mode 


difference 


0.106 


2.01537 2.01533 


0.00004 


0.312 


2.11152 2.01139 


0.00013 


0.507 


2.00710 2.00688 


0.00022 



the trapped BEC system. The merit of the trapped BEC experiment is that we can 
change the magnitude of the couphng constant by using the Feshbach resonance.^) 
Thus, the effect of the zero-mode may be experimentaUy observed with other set of 
parameter values. In this section, we investigate the couphng constant dependence 
of the zero-mode. 

We again assume an one-dimensional system and consider a BEC of rubidium 
atoms, ^'^Rb. The parameter values are the same as in S ^4.1I Here, we simply 
calculate the first excitation energy, which corresponds to the dipole mode. This 
is an ideal one-dimensional calculation, and it does not correspond directly to the 
experiments. 



Table II. Temperature and coupling constant 
dependences of the first excitation energy 
including the zero-mode. 





first excitation 


energy 


T/T, 


g^O.Ol 


g=0.015 


g=0.02 


1.5 X 10"^ 


1.0033 


1.0454 


1.0749 


3.7 X 10^2 


1.0685 


1.1243 


1.1908 


7.3 X 10-2 


1.143 


1.243 


1.350 



Table III. Temperature and coupling constant 
dependences of the subtraction of the fre- 
quency without the zero-mode from that 
that with the zero-mode. 



Table ^ lists the temperature and 
coupling constant dependences of the 
first excitation energy including the zero- 
mode. We can see that this energy be- 
comes larger as we increase not only the 
temperature but also the coupling con- 
stant. The difference*^ is listedd in Ta- 
ble Hn **) However, our numerical re- 
sults exhibit deviation from Kohn's theo- 
rem, according to which the energy of 
the dipole mode should be 1. This seems 
to be due to the fact that the dynamics of the non-condensate field are not prop- 
erly taken into account. As a matter of fact, the deviation becomes larger as the 
temperature is increased. Thus, the effect of the zero-mode observed here is not 







difference 






g^O.Ol 


g=0.015 


g=0.02 


1.5 X 10^2 


0.0007 


0.0020 


0.0037 


3.7 X 10-2 


0.0022 


0.0046 


0.0088 


7.3 X 10-2 


0.004 


0.008 


0.013 



*-* The definition of this quantity is the same as that given in Table Q 

**' The coupling constant g used here is smaller than 0.03, the value used in the calculations in 
the preceding subsections. Note that the spatial dimension is 1 in this subsection, while it is 3 in 
S^O audio 
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Table IV. Temperature and coupling constant 
dependences of the second excitation en- 
ergy including the zero-mode. 





second excitation 


energy 




g^O.Ol 


g=0.015 


g=0.02 


1.5 X 10"^ 


1.7874 


1.7873 


1.7974 


3.7 X 10-2 


1.8199 


1.8409 


1.8744 


7.3 X 10-2 


1.875 


1.928 


1.991 



Table V. Temperature and coupling constant 
dependences of the difference (second exci- 
tation energy). 



As we have discussed, the first ex- 
citation energy should be 1 for all tem- 
peratures and interactions, as asserted by 
Kohn's theorem. In order to determine 
the effect of the zero-mode, we should cal- 
culate the energy shift of the second exci- 
tation energy. Table HVI lists the temper- 
ature and coupling constant dependence 
of the second excitation energy including 
the zero-mode, and the difference is listed 
in Table IVl It is seen that the effect of the zero- mode is enhanced as the temperature 
and the coupling constant are increased. However, by comparing with Table 11111 we 
see that the effect of the zero-mode on the second excitation energy is of almost the 
same order as that on the first excitation energy. Thus, there is a possibility that 
this effect vanishes when we make corrections in order to satisfy Kohn's theorem. In 
short, the effect of the zero- mode in the one-dimensional system is not yet clear. 







difference 




T/T, 


g=0.01 


g=0.015 


g=0.02 


1.5 X 10-2 


0.0007 


0.0018 


0.0030 


3.7 X 10-2 


0.002 


0.004 


0.008 


7.3 X 10-2 


0.004 


0.007 


0.012 



§5. Concluding remarks 



We have formulated the dynamical response of a trapped BEC including the 
zero-mode at low temperature, using the HFBP and linear approximations. The 
zero-mode appears as a NG mode associated with the global phase symmetry break- 
ing. To take the zero-mode into account, we introduced the generalized Bogoliubov 
transformation following Ref. 20). We numerically calculated the dynamical response 
of the condensate by applying a time-dependent perturbation in the linear response 
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approximation. Thus, the fluctuations from the mean-field approximation are par- 
tiahy taken into account. The frequency of the condensate is determined by the 
first excitation energy of the quasiparticle. We observed that the presence of the 
zero-mode increases the frequency and hence the first excitation energy. With the 
parameter values used in the experiment/^) the effect of the zero- mode on the first 
excitation energy is very small. Thus, we conclude that the properties of the trapped 
BEC in this region can be described by the Boboliubov prescription, although it leads 
to an inconsistent description in quantum field theory. 

To determine the temperature and coupling constant dependences, we studied 
the first and second excitation energies in a one-dimensional system. Prom Kohn's 
theorem, it can be proved that the first excitation energy should be 1 for any temper- 
ature and interaction. Our results, however, exhibit deviation from 1. This deviation 
appears to be due to the inadequate approximation used for the non-condensate field, 
because the deviation becomes larger as the temperature is increased. The second 
excitation mode is also affected by the zero-mode. The magnitude of the energy shift 
is, however, almost the same as that of the first excitation energy. Thus, it is not 
clear whether the observed effect comes from the effect of the zero-mode or from the 
inadequate approximation violating Kohn's theorem. 

To improve our results, we should take account of the fluctuations of the non- 
condensate field, which is dropped in Eq. H3-9() . As a matter of fact, it has been 
reported that the collective dynamics of the non-condensate field play an important 
role.^^) In our calculation, if we take the fluctuations of the non-condensate fleld 
into account, the approximation is equivalent to the time-dependent linear Hartree 
approximation, that is, the RPA, and we can expect that a collective mode like the 
plasmon in an electron gas will be included. However, the translational invariance 
is broken by the trapped potential in this system. Thus, the analysis of the collec- 
tive mode is not simple, and we must introduce several approximations to simplify 
calculations.^^)' '^^^ The study of the effect of the fluctuations of the non-condensate 
fleld is a future problem. 

In this paper, we have concentrated on elucidating the effect of the zero-mode 
on the flrst excitation energy. However, calculations of the temperature dependence 
of the condensate fractions and the speciflc heat are also interesting. 

The effect of the gapless mode was also studied in Ref. 6). However, that work 
does not include the zero-mode in the non-condensate fleld. Consequently, the non- 
condensate fleld has zero projection on the condensate, and this property simplifies 
the calculation. By contrast, the zero-mode is included in the non-condensate field 
in our formalism. Thus, the gapless mode investigated in Ref. 6) is different from 
the zero-mode investigated here. 
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Appendix A 

Generalized Bogoliubov Transformation Including the Zero-Mode 

We have introduced the orthonormal set {wn{x)} satisfying Eq. (|2-15|) with the 
orthonormal and completeness conditions (|2-16a() and H2-16b|) . It is easily seen that 



Eq. (|2T5|) reduces to the GP equation for n = in the limt e ^ 0, because we have 
eo = 0. The condensate field is given by 

v{x) = yW^woix). (A-1) 

Using the orthonormal set, we expand (p{x) and introduce the creation and annihi- 
lation operators: 

oo 

ip{x) = ^anWn{x). (A-2) 

n=0 

The operators a„ and ali are subject to the condition 

[an,al/] = ^nn', (A-3) 

with all other commutators vanishing: [a„,a„'] = [an,aj^,] = 0. 

Substituting Eq. HA-2() into the mean field Hamiltonian (|2-13|1 with the condition 
(EUD), we obtain 

oo oo 
-ffMFA = ^(Cn + £e)aian + ^ ^2Unn'aian' + Unn'O-nan' + Unn'O^a^, ,(A-4) 
n=0 n,n'=0 

where 

f^in' = I j d^XWn{x)Wn'{x)v^{x). (A-5) 

The physical Fock space is defined by the creation and annihilation operators 
that diagonalize the mean field Hamiltonian. For this purpose, the following opera- 
tors are introduced: 

= J^^^-^(a. + at), (A-6a) 
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Pn = -m — 7y — K - a;)- 



These operators satisfy the CCRs 

[Qn,Pn] = iSnn'- 

All other commutators vanish. 

Then, Hmfa can be rewritten as 

n=0 n,n'=0 n=0 



(A-6b) 



(A-7) 



+ ee) + Un 



(A-8) 



Here, the matrix W has a block structure 
W = 



4(ee)C/oo + O(e^) V^^u'^ + 0(e 

V£^u' + 0{£^) W' + 0{£) 



(A-9) 



with 



u = 



( 4^C/io \ 



V 



(A- 10a) 
(A-lOb) 



(n, m = 1, 2, • 



This symmetric matrix W can be diagonalized by using an orthogonal matrix O with 
eigenvalues E"^ (n = 0, 1, 2, • • • ): 



with 



OWO'^ = E'^ 



{E)nm — E^djim • 



(A-11) 



(A-12) 



The above expressions are expanded with respect to the infinitesimal parameter 
£. It should be noted that the zero-th eigenvalue is given by 



(A-13a) 
(A- 13b) 



EQ = ^/ee\/EQ + 0{e^), 
Eo = AUoo - u^w'-^v! , 



and it vanishes when we take the limit e — >^ 0, but we keep it finite in the quantization 
procedure. 

Then, the orthogonal matrix O can be written as 



/ilO'W^'-V + 0(£2) 



O' + 0(£) 



,(A-14) 
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where O' is another orthogonal matrix diagonahzing the matrix W', 
with the diagonal matrix 



E 



/2 



El5nm + 0{e^). (n,m= 1,2,---) 



(A-15) 



(A-16) 



Thus the mixing among d„ (n = 1,2, • • • ) is regular with respect to s ; only the 
mixing between ao and other a„ gives rise to a singularity. 

By using the orthogonal matrix O, we can introduce new pairs of canonical 
operators as 



Qn — ^ ^ ^nmQn 



(A- 17a) 



m=0 



nmPm j 



m=0 



with [Qn,Pn'] = i^nn'j while all other commutators vanish. 
Finally, we can diagonalize the mean field Hamiltonian, 

°° 1 1 1 

i^MFA = Yl [2^n + ^Ql -2^^n + Se) + C/„ 
n=0 

oo ^ ^ 

^[Enbibn + -En ~ -(e^ 



ee) + Un 



n=0 



(A- 17b) 



(A-18) 



In the second equality, the operators Qn and P„ are related to 6„ and bn as 



Qn = J^ibn + bi), 



(A- 19a) 



Pn = ~t\[^{bn + hi) 



Thus, the physical vacuum is defined as 

S„|o) = 0. 



(A- 19b) 



(A-20) 



In short, the diagonalization studied to this point is realized by introducing the 
generalized Bogoliubov transformation. 



oo 

bn ~ ^ ^ j^C'nm'^m ~l" Snm^ 



m=0 

oo 



b\i — ^ ^ j^Cnrra^m "I" Snm.O'm 



(A-21a) 
(A-21b) 



m=0 
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Here the real matrices C and S in the matrix notation are given by 

C =- (E^Oe-^ + ^-^Oel) , (A-22a) 



2 

5 = i [E^Oe-^2 - E-^Oe^^ , (A-22b) 
with E defined by ()A-12|) and 

(e)nm = (en + ^^)^nm ■ (A-23a) 

One can easily check the properties of C and S from Eq. HA-22() . and we find 

oo 

^ ^ {CjiYnC^i'm ^nm^n'm) ^nn' ; (A'24a) 
) 

oo 

^ ^ (CnmSn'm ~ SnmCn'm) — • (A-24b) 



m=0 

oo 



m=0 



Thus, the field operators can be expanded in terms of 6-operators diagonalizing the 
Hamiltonian -f^MFA as 

oo 

^{x) = Y,[bne-'''-'wcn{x) - bie^'-'wsnix)] , (A-25a) 

n=0 
oo 

(^t(x) = Y,Pne'^"'wcn{x) - h.e-'^"'wsn{x)] , (A.25b) 



n=0 



with 



WCnix) = ^ CnmWm{x) , (A-26a) 

m=0 
oo 

■WSn{x) = ^ SnmWmix) . (A-26b) 
m=0 

It is emphasized that the canonical commutation relation 

[ip{x,t),^p\x',t)] =6^^\x-x') (A-27) 

holds in this quantization scheme. This is because we introduced the artificial inter- 
action l\'2-'6\i so that we could control the infrared divergence, as seen in Eqs. (|A-13|) . 
(fCTjl and 
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